Aspects of short range correlations in a relativistic model 
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In the present work short range correlations are introduced for the first time in a relativistic 
approach to the equation of state of the infinite nuclear matter in the framework of the Hartree-Fock 
approximation using an effective Hamiltonian derived from the a — lo Walecka model. The unitary 
correlation method is used to introduce short range correlations. The effect of the correlations in 
the ground state properties of the nuclear matter is discussed. 

: 

Nuclear physics is an effective theory of the nucleus regarded as a system of nucleons. In this theory the essential 
degrees of freedom are the center of mass coordinates of the nucleons as well as their spins and isospins, the interactions 
being expressed as nucleon-nucleon forces. It is usually accepted, however, that the nucleon-nucleon interaction 
(3JT)' becomes strongly repulsive at short distances in the relative coordinate of two nuclear particles. It is seen that 
the phase shifts for 1 So and 3 Si are positive at low and becomes negative at higher energies This indicates a 
' repulsive core at short distances and attraction at long distances. There are attempts to derive the nucleon-nucleon 
force using chiral perturbation theory 0]. In this approach, one- and two-pion exchange contributions are taken into 
account up to the third chiral order. However, in order to reproduce the nucleon scattering data and the D-wave 
phase shifts by this method, an ad hoc contact interaction (which represents the short range force) must also be 
included. The realistic nucleon-nucleon forces are basically phenomenological. The Bonn interactions Q are based 
^ ' on meson exchange which is treated in a relativistic non-local manner. The Argonne interactions Q, on the other 
hand, describe the pion exchange in a local approximation, the short and medium range nuclear interactions being 
' controlled by phenomenological parameters. 

Nonrelativistic calculations based on realistic NN potentials predict equilibrium points that are not able to describe 
simultaneously the correct binding energy and saturation density; either the saturation density is correct but the 
binding energy is too small, or the correct binding energy is obtained at a too high density. This behaviour is 
normally referred to "Coester line" This problem is generally circunvented through the introduction of a three- 
body repulsive force 0, or density-dependent repulsive mechanisms. A mechanism of this kind is already present 
in relativistic models. Due to Lorentz covariance and self-consistency, as the nuclear density increases, the nucleon 
effective mass decreases. As a result there is a reduction of the attractive force and a net increase of the repulsive 
force. In fact, the relativistic mean field (RMF) theory formulated by Teller and the others and by Walecka Q is 
successful in describing both infinite nuclear matter and finite nuclei. Moreover, the relativistic Hartree-Fock (HF) 
approximation [^.llOfl has also been used for the description of finite nuclei and infinite nuclear matter. 

In this paper, we have focused on the treatment of short range correlations in dense nuclear matter for the first 
time using relativistic equations of motion. Although there are several procedures which may be used to introduce 
short range correlations into the model wave function, we have preferred to work with the unitary operator method 
as proposed by Villars |ll|. There are several advantages in using an unitary model operator. In particular, one auto- 
matically guarantees that the correlated state is normalized. Further, we shall see that for calculation of correlation 
corrections to the matrix elements of one- and two-body operators is a simple task. The general idea exists for a long 
time 0,0,0] but has not been pursued for the relativistic case. This has been so because in non-relativistic models 
the interaction arises from the interplay between a long range attraction and a very strong short range repulsion so 
that, indeed, it is indispensable to take short range correlations into account. In relativistic mean field models, the 
parameters are phenomenological, fitted to the saturation properties of nuclear matter. Short range correlation effects 
are clearly included in the model parameters. We expect that the values of these parameters are more 'fundamental', 
or less artificial, when correlations are not neglected. 
We start with the effective Hamiltonian as 

H = J ipl(x)(-ia-'V + /3M) a pipp(x) dx 

+ ]^ J i'l l {x)^\{v)V a f}, 1 s{\x- y\)ips(y)il>i3(x) dx dy (1) 
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where 

„2 ~m a r 2 p -m^r 

V„(r) = -T- . V u> (r) = T- > ( 3 ) 

Air r Ait r 

and a are the Dirac-matrices. In the above, ip is the nucleon field interacting through the scalar and vector potentials. 
The equal time quantization condition for the nucleons is given by 

[Tl) a (£,t),ipp(y,tf] + =S a p5{x-y), (4) 

where a and j3 refer to the spin indices. We now also have the field expansion for the nucleons ip at time t=0 given 
as 151 



\= £ [U r {k)c r % + V r {-k)c\_ J e** (5) 
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For free spinor fields, we have cosx(/c) = M/e(k), sinx(fc) = \k\/e(k) with e{k) = \/k 2 + M 2 . However, we will 
deal with interacting fields so that we take the ansatz cosx(fc) = M*(k)/e*(k), sinx(fc) = \k*\/e*(k), with e*(k) = 



-•2 

k* + M* z (fc), where k* and M*(k) are the effective momentum and effective mass respectively. The equal time 

anti-commutation conditions are [c g, c = S rs Sr 2> — [c c, c -]+ . The vacuum | 0) is defined through 

c r | 0) = - | 0) = 0; one-particle states are written |fc, r) = c g | 0); two-particle and three-particle uncorrelated 

states are written, respectively as |fc, r; fc', r') = $ -. c, | 0), and \k, r; k 1 , r'; k", r") — -. c\ - c\ r - | 0), and so 
on. 

We now introduce the short range correlation through an unitary operator method. The correlated wave function 
is ^(J l^ 7 ) = e lS |$) where |$) is a Slater determinant and 5 is, in general, a n-body Hermitian operator, splitting 
into a 2-body part, a 3-body part, etc.. Consider now the expectation value of H, 

= (gjgjg) ($|e~^ g e ' g |j>) 

(*|*) ($|$) ' [ ' 

At this point we restrict ourselves to the two-body correlation diagrams shown in figure^ Let us denote the two-body 
correlated wave function by 



\k, r- k', r') = e lS \k, r; k', r') « f 12 \k, r; k' , r') (8) 

where /12 is the short range correlation factor, the so-called Jastrow factor |T^ . For simplicity, we have written 
/12 = /(r*i2), r*i2 = ri — r*2. We take the assumption that f(r) = 1 — (a + /3r) e _7r where a, /3 and 7 are parameters. 

At this point some remarks on the choice of the Jastrow factor are appropriate. The important effect of the short 
range correlations is the replacement, in the expression for the ground-state energy, of the interaction matrix element 

(k, r; k' , r'\Vi2\k, r; k' , r') by (k, r; k' , r'\V\2 + t\ + t2\k, r; k'.r') — (k, r; k', r'\U + t2\k, r; k' , r'), where U is the kinetic 
energy operator of particle i. As argued by Moszkowski [ijj and Bethe it is expected that the true ground-state 
wave function of the nucleus containing correlations coincide with the independent particle, or Hartree-Fock wave 
function, for interparticle distances r > rheal > where Theai ~ 1 fm is the so-called "healing distance" . This behavior is 
a consequence of the constraints imposed by the Pauli Principle. Moreover, although in general the correlation factor 
f(r) may depend on the isospin and spin quantum numbers of the two-body channel, we assume, for simplicity, that 
it is a plain, state independent Jastrow factor. A natural consequence of having the correlations introduced by an 
unitary operator is the occurrence of a normalization constraint on /(r), 

(/ 2 (r) - 1) d 3 r = 0. (9) 
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The correlated ground state energy becomes 

£ = 4 r k 2 dk [\k\sm X (k) + Mcos X (fc)] + EMfi + ^Mfe 
^ Jo 2 2 
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+ Cops^ I k 2 dk [\k\sinx{k) + Mcos X (k)} 
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' k 2 dkk' 2 dk' | |fc|sinx(fc) + 2 Mcosx(fc)J/(fc,fc') + |fc| sinx(fc') J(fc,fc')} 

Ai(fc,fc') +cosx(fc)cosx(fc') ^ B^fc,^) +sinx(fc)sinx(fc') ^(M') 



fc dfc fc' d/c' 



(10) 



where, Co — J (f 2 {r) — 1) d 3 r so that, according to @, Co = 0, and Aj, £?i, (7, , / and J are exchange integrals. 
In the above equation for the energy density, the first term results from the kinetic contribution, the second and 
third terms come respectively from the a and u> direct contributions from the potential energy with correlations, 
the fourth and fifth ones from the direct and exchange correlation contribution from the kinetic energy, and the last 
one from the a + u) exchange contributions from the potential energy with correlations. The angular integrals are 
given by Ai(k,k') = B t {k,k') = 2tt gf/An £ dcos 6 F t {k, k' , cos 9), C t {k,k') = 2tt gf/Air Jfi cos 6 d cos 6 F t (k, k', cos 0), 
I(k, k') = 2tt Jq dcos9 C\{k, k', cos 9), and J(k, k') = 2ir J™ cos 9 dcos9 C\{k, k', cos 9), where 

Fi(k,k') = I [f(r)V T {r)f(r)] dr and C r (k, k') = / (/ 2 (r) - 1) e ^- n '> r dr. (11) 



More explicitly, the first terms of the above angle integrals read for f(r) = 1 

^cr = gl^cr B a = g 2 9 a C a = -2g 2 (\> a 



= 2g 2 J UJ B u) = -±g 2 J u) C^ = -4g 2 



. ,' (fc + fcQ 2 +mf \ , k 2 + k' 2 + mi 2 Q 



The baryon density and the scalar density are 



2 v f k f - 2 v k 3 f 2v f k f - - 

PB = TTTTT / dk = a ? ■ />« = TTTTT / cosx(fc) dfc. (12) 
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The parameters of the model have to be fixed. They are the couplings g a , g Ul the meson masses, m a and m u and 
also three more parameters from the short range correlation function, a, f3 and 7. The couplings are chosen so as to 
satisfy the ground state properties of the nuclear matter and are given in table Q] We choose m a — 550 MeV and take 
th u = 783 MeV. The normalization condition © determines (3. We fix a either by imposing the condition /(0) = 0, 
which appears to be a natural choice from our experience with the non-relativistic case, or by minimizing the energy. 
We choose 7 so that the correct healing distance 0] is reproduced. If a is chosen to be 1, we assume a density 
independent parameter 7 equal to 750 (HF+corr-IV) . On the other hand, if we choose to determine a variationally, 
we assume that the parameter 7 increases linearly with the Fermi momentum. The last choice is consistent with 
the idea that the healing distance decreases as kp increases. Of course there are other possible choices for these 
parameters. The parameters we have used are tabulated in table [I] together with the compressibility K, the relative 
effective mass Ad* /M, the kinetic energy T/ps — M, the direct and exchange parts of the potential energy (Vd/pB 
and V e /pB repectively) with correlation and the correlation contribution to the kinetic energy T c / P b 1 all computed 
at the saturation point. We next compute the range of the correlations, obviously related with the healing distance, 
also included in tableland defined as [l9j : 

p2 _ / r 2 {f{r)-l) 2 d*r 

f(f(r) - l) 2 d*r ■ [ 6 > 

In figure HI we plot the binding energies as function of the density for the Hartree, HF, quark-meson-coupling model 
(QMC) [20|. a non-linear Walecka model NL3 and the 4 choices of the parameters in our calculation, as given in 
table [I] From this figure one can see that in all cases that inclusion of correlations make the equation of state (EOS) 
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softer than Hartree or HF calculations. NL3 and QMC also provide softer EOS around nuclear matter saturation 
density but around two times saturation density, some of the EOS with correlations are softer than NL3. Correlations 
always tend to soften the EOS, except when a — 1 is fixed. In this case, the EOS coincides with the curve obtained 
with the Hartree calculation. 

In figure |3] we plot the individual contributions to the energy density in function of the density for the Hartree, 
HF and HF+correlations (set III). Notice that the correlation contribution is of the same order as the exchange 
contribution. Hence, it cannot be disregarded. In fact if we compare the coupling constants g a and obtained 
for the different calculations we conclude that the introduction of the correlations reduces the coupling constants. 
Correlation effects in the Hartree and HF calculations are taken into account by a correct choice of the coupling 
constants. However, the explicit introduction of correlations has other effects such as softning the EOS. 

In conclusion, we have included two-body correlations in the a — uj relativistic model. We have seen that they 
correspond to an important amount of the energy density and provide a softer EOS as compared with the HF 
approximation. Furthermore, the effect of the correlations is similar to the role played the non-linear terms in 
relativistic models or the 3-body forces in non-relativistic potential models. The inclusion of three-body correlations 
is currently under investigation. 

This work was partially supported by CNPq (Brazil), CAPES (Brazil), GRICES (Portugal) under project 100/03, 
FEDER and FCT (Portugal) under the projects POCTI/FP/FNU/50326/2003 and POCTI/FIS/451/94. 



[1] M.H. MacGregor, A. R. Arndt, and R.M. Wright, Phys. Rev. 182, 1714 (1969) and references therein. 
[2] D.R. Entem and R. Machleidt, Phys. Lett B 524, 93 (2001). 

[3] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989); R. Machleidt, Phys. Rev. C 63, 024001 (2001). 
[4] R. Wiringa, R. Smith, and T. Ainsworth, Phys. Rev C 29, 1207 (1984). 

[5] F. Coester, S. Cohen, B.D. Day, and CM. Vincent, Phys. Rev. C 1, 769 (1970); R. Brockmann and R. Machleidt, Phys. 
Rev. C 42, 1965 (1990). 

[6] R.B. Wiringa, V. Fiks and A. Fabrocini, Phys. Rev. C 38, 1010 (1988); W. Zuo, A. Lejeune, U. Lombardo and J.-F. 

Mathiot, Nucl. Phys. A 706, 418 (2002). 
[7] H.P Diirr and E. Teller, Phys. Rev. 101, 494 (1956); M. H. Johnson and E. Teller, Phys. Rev. 98, 783 (1955); H.P Diirr, 

Phys. Rev. 103, 469 (1956). 

[8] B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986); J.D. Walecka, Ann. of Phys. 83, 491 (1974); B.D. Serot, J.D. 

Walecka, Int. J. Mod. Phys. E6, 515 (1997). 
[9] C.J. Horowitz and B.D. Serot, Nucl. Phys. A399, 529 (1983). 
[10] A. Bouyssy, J.-F. Mathiot and N.V. Giai and S. Marcos, Phys. Rev. C 36, 380 (1987). 

[11] F. Villars, "Proceedings of the International School of Physics, 'Enrico Fermi'-Course 23, (1961)." Academic Press, New 
York, 1963; J.S. Bell, " Lectures on the Many-Body Problem, First Bergen International School of Physics." Benjamin, 
New York, (1962); F. Coester and H. Kiimmel, Nucl. Phys. 17 477 (1960). 

[12] J. da Providencia and CM. Shakin, Ann. Phys.(NY) 30, 95 (1964). 

[13] K. Suziki, R. Okamoto and H. Kumagai, Phys. Rep. 242 181 (1994). 

[14] H. Feldmeier, T. Neff, R. Roth, and J. Schnack, Nucl. Phys. A 632, 61 (1998); T. Neff and H. Feldmeier, Nucl. Phys. 713 
311 (2003). 

[15] A. Mishra, P.K. Panda, S. Schramm, J. Reinhardt and W. Greiner, Phys. Rev. C 56, 1380 (1997). 

[16] J. da Providencia and C. M. Shakin, Phys. Rev C 4, 1560 (1971); C. M. Shakin, Phys. Rev C 4, 684 (1971). 

[17] R. Jastrow, Phys. Rev. 98, 1479 (1955). 

[18] S.A. Moszkowski and B.L. Scott, Ann. Phys. (N.Y.), 11, 65 (1960). 
[19] H. Bethe, Ann. Rev. Nucl. Sci. 21, 93 (1971). 

[20] P. A. M. Guichon, Phys. Lett. B 200, 235 (1988). K. Saito and A.W. Thomas, Phys. Lett. B 327, 9 (1994); P.K. Panda, 

A. Mishra, J.M. Eisenberg, W. Greiner, Phys. Rev. C 56, 3134 (1997). 
[21] G. A. Lalazissis, J. Konig and P. Ring, Phys. Rev. C 55, 540 (1997). 



■5 



TABLE I: Model parameters and ground state properties of nuclear matter at saturation density are given. These results were 
obtained with fixed: M = 939 MeV m CT = 550 MeV, = 783 MeV at k F0 = 1.3 fin -1 with binding energy E B = e/p - M = 
— 15.75 MeV. We have used 4 choices of the parameters for the correlation. HF+corr-I corresponds to 600 + 200&f/&fo, 
HF+corr-II to 700 + 200k F /k F0 , HF+corr-III to 600 + 300fc F /fc F0 and HF+corr-IV corresponds to a = 1 and 7 fixed. 
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FIG. 1: The kinetic and the potential contribution to the two-body cluster diagrams 




FIG. 2: EoS for different parametrizations as defined in the table. 
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FIG. 3: Individual contributions for the energy density 



